Effective pseudopotential for energy density functionals with higher order derivatives 
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We derive a zero-range pseudopotential that includes all possible terms up to sixth order in 
derivatives. Within the Hartree-Fock approximation, it gives the average energy that corresponds 
to a quasi-local nuclear Energy Density Functional (EDF) built of derivatives of the one-body density 
matrix up to sixth order. The direct reference of the EDF to the pseudopotential acts as a constraint 
that divides the number of independent coupling constants of the EDF by two. This allows, e.g., 
for expressing the isovector part of the functional in terms of the isoscalar part, or vice versa. We 
also derive the analogous set of constraints for the coupling constants of the EDF that is restricted 
by spherical, space-inversion, and time-reversal symmetries. 
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I. INTRODUCTION 

One of the big challenges of the current research in nu- 
clear structure physics is the search for a universal energy 
density functional (EDF) Among different possible 
approaches to this search, the consideration of a local or 
quasi-local EDF based on the density-matrix expansion 
QDME) is in recent years the object of intense studies 
[^-lZ] ■ These aim at improving the classic work of Negele 
and Vautherin and better theoretical understanding 
based on the effective theory [HI and on the frame- 
work of the density functional theory [l^ . 

In the recent work , we proposed a new expansion of 
the nuclear energy density in higher-order derivatives of 
densities. There, following the effective-theory approach, 
a Skyrme-like quasi-local next-to-next-to-next-to-leading 
order (N^LO) EDF was derived with terms of the EDF 
constrained only by symmetry principles. In the present 
derivation, we took the route in opposite direction as 
compared to what has been done for the standard Skyrme 
next-to-leading order (NLO) EDF. Namely, historically, 
the Skyrme force has been initially proposed first as an 
expansion of the effective interaction in relative momenta 
up to second order [U, Next, for this force the aver- 
age Hartree-Fock (HE) energy was evaluated, giving the 
Skyrme EDF with half of the coupling constants con- 
straint to the other half, see Ref. [l5| for the modern 
complete analysis. Only later, a possibility of releasing 
these constraints was considered and studied, see, e.g., 
Ref. [l^ for the analysis of the spin-orbit term. 

In the present work we complete the results of Ref. 
by deriving the expansion of the effective interaction 
in relative momenta up to N^LO. This generalizes the 
Skyrme force up to sixth order and allows us to make 
a link with the general N-^LO EDF derived in Q. One 
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should stress that the present analysis is not at all an 
independent repetitive derivation of the same functional. 
Indeed, the constraints on the EDF coupling constants, 
which are induced by the HE averaging of this general- 
ized force, cannot be obtained without following the path 
presented in this study. 

The complete higher-order EDFs or pseudopotentials 
have never yet been applied in practical calculations. The 
work towards this goal is now in progress, with basic 
derivations like the ones of Ref. and in the present 
work coming first, the construction of numerical codes 
like the one in Ref. [17] coming next, and the full ad- 
justments of coupling constants that will follow. In this 
respect, at present, we are in a similar phase of stud- 
ies as before the chiral N'^LO potentials for two-nucleon 
systems were adjusted, see, e.g., Ref. [l8l|, and after the 
tools for calculating the corresponding N'^LO diagrams 
were developed, see, e.g., Ref. [l^. Nevertheless, studies 
of particular higher-order EDF terms have already been 
performed 0, |2l| . 

In Ref. [13], the question of convergence of the series 
in higher-order derivatives was recently addressed within 
the DME applied to the Gogny non-local functional, and 
it was shown that every next order up to sixth gives con- 
tributions smaller by large factors. This gives us confi- 
dence that fits of higher-order EDFs have a fair chance 
of converging. A rigorous power counting scheme, anal- 
ogous to what has been introduced in the chiral pertur- 
bation theory [2^, would have to use derivatives of reg- 
ularized zero-range interactions, see, e.g. Ref. [l3|- Such 
a regularization would provide a proper cut-off scale, 
against which the powers of derivatives could be esti- 
mated. A good model of the regularized delta force is 
the Gaussian interaction, which, however, leads (through 
the exchange term) to non-local functionals. Within the 
EDF methodology, an effective theory based on deriva- 
tives of finite-range force is in principle possible, but has 
not yet been tried because of the degree of numerical com- 
plications involved. In the language of the effective field 
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theory, the power counting scheme allows us to properly 
classify diagrams of the perturbation series, however, the 
ideas of an effective theory are much more general than 
their applications in the field theory - here we use them 
within the framework of standard quantum mechanics of 
many-body systems. 

The EDF description of nuclear states is phenomeno- 
logical in the sense that it depends on the coupling con- 
stants, which are usually fitted to available experimen- 
tal data, see recent Refs. [13, HEl on fitting the second- 
order (NLO) Skyrme functionals. Fits of the full N'^LO 
EDF are much more complicated because of strong inter- 
dependencies of the coupling constants and instabili- 
ties 26] occurring in certain regions of the parameter 
space. Our main motivation to carry out the present 
work was to find constraints on parameters of the general 
EDF, which result from its relation to a pseudopotential. 
Such a relation reduces the number of parameters that 
have to be fit to data, and by this virtue is a positive 
change, at least at the preliminary stage of adjustments. 

Instead of fitting the coupling constants of the EDF, 
it is also possible to derive them directly using the 
DME [H. The DME gives an EDF which approxi- 
mates more complicated and time consuming HF cal- 
culations based on finite-range forces. When applying 
the DME, the relations to pseudopotentials are however 
usually broken Q. By enforcing these relations, as done 
here, one ensures that the generated EDF is free from 
unphysical self-interaction [23, and can be applied 
in beyond-mean-field applications without problems, see, 
e.g., Refs. (Hill. 

By following the standard convention, here we call the 
generalized Skyrme force pseudopotential, which is the 
name denoting a quasi-local operator depending on spa- 
tial derivatives. We also consequently use the names 
'parameters' to denote numerical coefficients of different 
terms of the pseudopotential, and we use the names 'cou- 
pling constants' to denote numerical coefficients of terms 
in the EDF. 

The paper is organized as follows. In Sec. |lT] we con- 
struct the pseudopotential in two alternative forms and 
list all its terms up to N^LO. We also evaluate the con- 
straints imposed by the gauge symmetry. In Sec. IIIII we 
discuss the procedure of HE averaging to obtain the EDF 
from the pseudopotential. In particular, in Scc. lIII Al we 
derive the general relations connecting the parameters 
of the Galilean-invariant pseudopotential to the coupling 
constants of the EDF, whereas in Sec. lIIIBl we derive the 
constraints for the case of conserved gauge symmetry. In 
Sec. lIVI we reduce our results to the case of the conserved 
spherical, space-inversion, and time-reversal symmetries. 
After formulating the conclusions of the present study 
in Sec. El in Appendices El IC] we present derivations 
related to the time-reversal invariance and hermiticity 
of the pseudopotential, we list results pertaining to the 
gauge-invariant pseudopotentials, and we give relations 
between the two alternative forms of pseudopotentials. 
Results obtained in the present work that are too volu- 



minous to be published in the printed form are collected 
in the supplemental material [3l|. 



II. GENERAL FORM OF THE 
PSEUDOPOTENTIAL IN THE 
SPHERICAL-TENSOR FORMALISM 

A. Central-like form of the pseudopotential 

The Skyrme interaction is one of the most impor- 
tant phenomenological effective interaction used in mi- 
croscopic nuclear structure calculations: such two-body 
interaction is a short-range expansion up to the second 
order in derivatives, which contains a certain number of 
fit parameters adjusted to reproduce the experimental 
data. In the literature the Skyrme interaction is usually 
written in cartesian representation, but for our extended 
pseudopotential we adopt the spherical-tensor represen- 
tation of operators [s^l, whose building blocks can be 
found in '5]. 

Depending on the specific form of the coupling of the 
derivative operators with the spin operators, different 
ways to construct the pseudopotential are possible. A 
particular form of the pseudopotential, which we call 
central-like or LS-like, is constructed in the present Sec- 
tion. It is based on coupling together the derivative oper- 
ators and spin operators, which are then coupled to rota- 
tional scalars. An alternative form, called tensor-like or 
JJ-like, is presented in Section HlDI There, each deriva- 
tive operator is coupled with one spin operator, and then 
they are coupled together to rotational scalars. 

In the central-like form, the pseudopotential is a sum 
of terms, 



/ ^ nL,vi20 nL,vi2 



(1) 



n L , 



where the sum runs over the allowed indices of the tensors 
according to the symmetries discussed below. Each term 
in the sum is accompanied by the corresponding strength 
parameter C^^^^ g, and explicitly reads. 
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X (^l-P*^P-p-j5i2Kri;rir2). (2) 

In Eq. ([2]), are the spherical tensor derivatives of 

order h and rank L built of the spherical representations 
of the relative momenta k — (Vi — V2)/2i, 



fci,M={-i,o,i} - 'XTS 

^ZJ V3 ^ ! 



(3) 
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TABLE I: Derivative operators KnL up to N^LO as expressed 
through spherical tensor representation of relative momenta 
k defined in Eq. Q. 



No. 



tensor KnL 



order n rank L 



1 
2 
3 

4 
5 
6 
7 

8 
9 
10 
11 
12 
13 
14 
15 
16 



1 
k 

[kk]^ 

[k[kkU, 

l2 



J2 

■J2J314 



[fe[fc[fcfc]j 

[kkUkikkU, 
[kk]l 

[fcfe]Q[fefc]2 

[kk\^[k[k[kku^\^ 

[k[k[k[k[kk]^]^],],;^. 



up to sixth order they are listed in Table HI Similarly, 
operators K'_ - are built of the relative momenta k' — 

nL 

(V'l - V'2)/2i. 

The symmetrized two-body spin operators Sy-^^s are 
defined as, 



(4) 



(i) 

where vi2 = vi + V2 and a^/i are the spherical-tensor 
components of the rank-u Pauli matrices acting on spin 
coordinates of particles i = 1 or 2. They are expressed 
as 



a 



1, ^={-1,0,1} 



= 1, 



(5) 



V2 



P)} (6) 



through the spin unity matrix 1 and the standard Carte- 

(i) 

sian components of the Pauli matrices cr^'y^z- 
The Dirac delta function, 

S 12 {r[r 2, rir 2) = S{r[-ri)S{r2-r2)S{ri-r2) 

= S{r[-r2)S{r'2-ri)S{r2-ri). (7) 



ensures the locality and zero-range character of the pseu- 
dopotential. The action of derivatives Kf^i and 

on Si2{r'ir2,rir2) has to be understood in the stan- 
dard sense of derivatives of distributions. Whenever the 



pseudopotential ([T]) is inserted into integrals to calculate 
the two-body matrix elements, the integration by parts 
transfers the derivatives onto appropriate wave functions 
in the remaining parts of integrands. 

The exchange term is explicitly embedded in the pseu- 
dopotential through the operator 



pMpapr ^ (-l)«'-(l + V3 



.(1)^(2) 



-V3 



^(1)^(2) 
'1 ~1 



^(1)^2) 



,(8) 



where r}*'' are the standard spherical-tensor isospin Pauli 
matrices defined analogously as in Eq. The square 
brackets with superscripts and subscripts denote the cou- 
pling of spherical tensors in the isospin space and co- 
ordinate space, respectively. The above definitions and 
conventions exactly correspond to those introduced in 
Ref. 0. 

The zero range of the pseudopotential has an impor- 
tant bearing on the structure of terms in Eq. ([2]). Indeed, 
only for the zero-range force, the space-exchange (Majo- 
rana) operator can be replaced, in any individual 
term, by the phase (—1)" appearing in Eq. ([5]). More- 
over, apart from the isospin-exchange operator P"^, terms 
of the pseudopotential cannot then depend on isospin. 
This fact, effectively reduces by half the number of al- 
lowed terms of the pseudopotential, as compared to what 
would have been possible for a finite-range potential. 
This is at the origin of the numbers of allowed terms 
of the pseudopotential being equal one half of the num- 
bers of the allowed terms of the EDF, which we discuss 
below. 

The full antisymmetrization of the pseudopotential in- 
cludes the exchange operator in the isospin space; there- 
fore, in the following we consider the EDF with the 
isospin degree of freedom included, that is, we discuss 
both the isoscalar and isovector terms of the N'^LO 0], 
which allows us to fully incorporate the pro ton- neutron 
mixing at the level of the energy density [l5|. 

The general form of the pseudopotential and the al- 
lowed terms listed below refiect the fact that the funda- 
mental symmetries of the two-body interaction must be 
respected, see Appendix[Xl In particular, (i) all terms are 
scalar operators, that is, they are coupled to the total an- 
gular momentum 0, which ensures the rotational invari- 
ance, (ii) the total number of derivative operators must 
be even, namely, n + h' = 0,2,4,6, which ensures the 
time-reversal and parity invariances, (iii) the parameters 
C"-^ ^ of the pseudopotential must be real, to guar- 

nL,vi2S > o 

antee both the time-reversal invariance and hermiticity, 
and (iv) the invariance under exchange of the coordinates 
of particle 1 and 2 is respected by expression 
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TABLE II: Zero-order terms of the pseudopotential ([2]). 



No. 


n' 


L' 


h 


L 


"12 


S 


gauge 


1 




















Y 


2 














2 





Y 


TABLE 
terms. 


III: 


Same as 


in Table [TI] but for the second order 


No. 


n' 


I' 


n 


L 


Vl2 


s 


gauge 


I 


2 

















Y 


2 


2 











2 





Y 


3 


2 


2 








2 


2 


Y 


4 


1 


1 


1 


1 








Y 


5 


1 


1 


1 


1 


2 





Y 


6 


1 


1 


1 


1 


1 


1 


Y 


7 


1 


1 


1 


1 


2 


2 


Y 



TABLE IV: Same as in Tablellllbut for the fourth order terms. 



No. 


n' 


L' 


n 


L 


?^12 


5 


gauge 


1 


4 

















D 


2 


4 











2 





D 


3 


4 


2 








2 


2 


D 


4 


3 


1 


1 


1 








Y 


5 


3 


1 


1 


1 


2 





Y 


6 


3 


1 


1 


1 


1 


1 


N 


7 


3 


1 


1 


1 


2 


2 


D 


8 


3 


3 


1 


1 


2 


2 


I 


9 


2 





2 











D 


10 


2 





2 





2 





D 


11 


2 


2 


2 





2 


2 


D 


12 


2 


2 


2 


2 








I 


13 


2 


2 


2 


2 


2 





I 


14 


2 


2 


2 


2 


1 


1 


N 


15 


2 


2 


2 


2 


2 


2 


I 



B. Lists of terms of the pseudopotential V order by 

order 

In Tables UniVl are listed, respectively, all possible terms 
of the pseudopotential ([T]) in zero, second, fourth, and 
sixth order. In each order, the numbers of terms equal 2, 
7, 15, and 26, giving the total number of 50 terms up to 
N'^LO. We see that these numbers of terms are exactly 
equal to those corresponding to the EDF in each isospin 
channel with the Galilean invariance imposed, cf. Ta- 
ble VI of Ref. ,^]. One should note that each term of 
the pseudopotential © is Galilean-invariant by construc- 
tion, because it is built with relative-momentum opera- 
tors Kf^i; therefore, the pseudopotential is not changed 
by a transformation to a system moving with a constant 
velocity. When both isoscalar and isovector channels are 
considered in the EDF, the number of EDF terms be- 
comes in each order twice larger than the number of terms 
of the pseudopotential. 

This means that the EDF obtained by averaging the 
pseudopotential is constrained by as many conditions as 
there are terms in each isospin channel. One possible 
solution is than to find a one-to-one correspondence be- 
tween the EDF and the pseudopotential by relating the 
isoscalar part of the EDF to its isovector part, in a way 
that will be showed explicitly in the following Sections of 
this work. 

To make the connection between the pseudopotential 
and the standard form of the Skyrme interaction more 
transparent, we give here the relations of conversion be- 
tween the parameters of the zero- and second-order pseu- 
dopotential and those of the Skyrme interaction, see 



Ref. [T^ for the definitions used. They read. 



to — C'oO.OO + ^^^00,20' (9^) 

toxo - -^C%o, (9b) 

^1 = "y^^OaOO + 2^*00 20, (9c) 

h^i = -^^20 20, (9d) 

*2 = —^Ciim + 2C'n,20' (9e) 

^2X2 = ~lcll^o^ (9f) 

Wo = ^^nai' (9g) 

to = -^^i"22, m 

te = ~ /;r^00.22- (9i) 



In relations of Eqs. ©, parameters t^ and t^xs are 
missing: they are related to the terms of the Skyrme 
interaction depending on density, which have been in- 
troduced to mimic the effects of the three-body force in 
the phenomenological interaction and to get the satura- 
tion feature of the nuclear force. In the same way, the 
zero-order parameters Cg" qq and C00.20 of the pseudopo- 
tential, see Eqs. (Pa)) and ((9b| . should become density- 
dependent. 

In his effective nuclear potential, Skyrme also intro- 
duced [T^ one additional term of the fourth order, which 
he justified through the presence of considerable D- 
waves in the nucleon-nucleon interaction energies around 
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TABLE V: Same as in Table HIl but for the sixth order terms. 



No. 


n' 


L' 


n 


L 


"12 


S 


gauge 


1 


6 

















D 


2 


6 











2 





D 


3 


6 


2 








2 


2 


D 


4 


5 


1 


1 


1 








D 


5 


5 


1 


1 


1 


2 





D 


6 


5 


1 


1 


1 


1 


1 


N 


7 


5 


1 


1 


1 


2 


2 


D 


8 


5 


3 


1 


1 


2 


2 


I 


9 


4 





2 











D 


10 


4 





2 





2 





D 


11 


4 


2 


2 





2 


2 


D 


12 


4 





2 


2 


2 


2 


D 


13 


4 


2 


2 


2 








I 


14 


4 


2 


2 


2 


2 





I 


15 


4 


2 


2 


2 


1 


1 


N 


16 


4 


2 


2 


2 


2 


2 


D 


17 


4 


4 


2 


2 


2 


2 


I 


18 


3 


1 


3 


1 








D 


19 


3 


1 


3 


1 


2 





D 


20 


3 


1 


3 


1 


1 


1 


N 


91 


Q 
O 


1 

_L 


Q 
O 


1 

_L 


9 


9 




22 


3 


3 


3 


1 


2 


2 


D 


23 


3 


3 


3 


3 








I 


24 


3 


Q 
O 


Q 
O 


Q 
O 


9 





I 


25 


3 


3 


3 


3 


1 


1 


N 


26 


3 


3 


3 


3 


2 


2 


D 


100 MeV. Also 


in this 


case, 


we £ 


nve the 


relation between 



the corresponding parameter and the parameter of 
our full pseudopotential, 



- 2^20,20- (10) 



C. Gauge invariance of the pseudopotential 

Besides the Galilean invariance mentioned above, the 
standard Skyrme force has been also proved to be invari- 
ant with respect to a more general local gauge invariance, 
and to give rise to the energy density that is invariant un- 
der the same symmetry when specific relations between 
the coupling constants are set [s^, Is^] • 

The gauge transformation acts on a many-body wave 
function by multiplying it with a position-dependent 



phase factor, that is, 

l*')-exp^z^0(r,)j 1*), (11) 
and its action transferred onto the pseudopotential is. 

Apart from zero order, the terms of the pseudopoten- 
tial are not trivially invariant with respect to the trans- 
formation of the Eq. (|12l) and, in general, the transformed 
pseudopotential V' is different than the original pseu- 
dopotential V. To impose the gauge invariance on the 
pseudopotential, one has to derive a list of constraints 
among the parameters, which can be done using the con- 
dition 

[^in),V] + [^ir2),V]=0. (13) 

As expected, at second order, all the 7 terms of the pseu- 
dopotential listed in Table Hill fulfill condition (|13p . Then 
they all are the stand-alone gauge invariant terms of the 
pseudopotential, which in the last column of the Table 
is marked by the letter Y. On the other hand, at fourth 
order, only two of the terms of the pseudopotential listed 
in Table IIVI those that correspond to parameters Cfl qq 
and Cfl fulfill condition p^ . At sixth order, none of 
the terms are stand-alone gauge invariant. 

At fourth order, the gauge invariance forces seven pa- 
rameters of the pseudopotential to be specific linear com- 
binations of four independent ones. In Table ITVl they are 
marked by letters D and I, respectively. In Appendix IB| 
we list such relations between the dependent and inde- 
pendent parameters. One should note that other choices 
of the four independent parameters are also possible, that 
is, at fourth order, there are simply four different gauge- 
invariant linear combinations of terms of the pseudopo- 
tential ([Ij. Moreover, at this order, there are also two 
terms that alone are gauge non-invariant - those that 
correspond to parameters Cfl i^ and CH^; in Table HVl 
they are marked by letters N. Similarly, at sixth order, 
there are six gauge-invariant linear combinations of terms 
of the pseudopotential, that is, sixteen dependent param- 
eters are related to six independent ones, see Appendix IbI 
and there are also four alone gauge non-invariant terms 
corresponding to parameters Cfl C||;^j, C^l n, and 
r'33 

^33,11- 

A comparison between the numbers of terms of 
the Galilean-invariant pseudopotential and the gauge- 
invariant pseudopotential is plotted in Fig. [TJ Again we 
note that at each order, the numbers of gauge-invariant 
parameters (2 for the zero order, 7 for the second order, 
6 for the fourth order, and 6 for the sixth order) are ex- 
actly the same as the numbers of independent coupling 
constants of the EDF in each isospin channel with the 
gauge invariance imposed, cf. Table VI of Ref. ^. A gam, 
this observation will be crucial when we proceed to derive 
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I Galilean-invariant pseudopotential 
> gauge-invariant pseudopotential 
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Order in derivatives 

FIG. 1: (Color online) Number of terms of the pseudopoten- 
tial ([2}, plotted as a function of the order in derivatives. 



relations between the isoscalar and the isovector parts of 
the EDF, stemming from the gauge- invariant pseudopo- 
tential. We also remark that whereas the second-order 
spin-orbit term, corresponding to parameter CHn, is 
gauge invariant, all higher-order spin-orbit terms, corre- 
sponding to parameters C^-^^^ with ; 
the gauge symmetry. 



D. Tensor-like form of the pseudopotential 

In this Section, we present the tensor-like form of the 
pseudopotential, which is, in fact, a different form of cou- 
pling of the relative-momentum operators with the spin 
operators, just like in the tensor term of the standard 
Skyrme interaction. In this form, the pseudopotential of 
Eq. ([Ij is a sum of the following terms. 



V 



where 



nL,vi2J 



Z ^ 7iL,vi2J nL,vi2J^ 

n L' . 



(l 2 ^1-1 ,1^2 ) >^ 



(14) 



n' L' ^1 



fi'U ""i 



M Act At 



l-P'''P^Pn6i2ir[r'^;r,r2). (15) 



The lists of the zero-, second-, fourth-, and sixth-order 
terms VT-^ , of the pseudopotential are given, respec- 

nL,vi2J ^ o 7 i 

lively, in Tables IVIHIXl which are the analogues of Ta- 
bles HIHY] given in Section HTB] 



TABLE VI: Zero-order terms of the recoupled pseudopoten- 
tial (nni. 



No. 


n 


L' 


n 


L 


V12 J 


1 

















2 














2 


TABLE VII: Same 
terms. 


as in 


Table IVII but for the second-order 


No. 


n' 


L' 


n 


L 


V12 J 


1 


1 


1 


1 


1 


1 


2 


1 


1 


1 


1 


1 1 


3 


1 


1 


1 


1 


2 


4 


1 


1 


1 


1 


2 1 


5 


1 


1 


1 


1 


2 2 


6 


2 














7 


2 











2 1 


8 


2 


2 








2 1 



By means of the recoupling technique, it is possible to 
determine relations between the two different coupling 
schemes of the pseudopotential. This derivation, along 
with the relationships between the corresponding param- 
eters C?/" ^ and C?/" is presented in Appendix O 



TABLE VIII: Same as in Table but for the fourth-order 



terms. 


No. 


n 


L' 


n 


L 


V12 


J 


1 


2 





2 











2 


2 





2 





2 


1 


3 


2 


2 


2 


2 





2 


4 


2 


2 


2 


2 


1 


2 


5 


2 


2 


2 





2 


1 


6 


2 


2 


2 


2 


2 


1 


7 


2 


2 


2 


2 


2 


2 


8 


2 


2 


2 


2 


2 


3 


9 


3 


1 


1 


1 





1 


10 


3 


1 


1 


1 


1 


1 


11 


3 


1 


1 


1 


2 





12 


3 


1 


1 


1 


2 


1 


13 


3 


1 


1 


1 


2 


2 


14 


3 


3 


1 


1 


2 


2 


15 


4 

















16 


4 











2 


1 


17 


4 


2 








2 


1 
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TABLE IX: Same as in TablelVTlbut for the sixth-order terms. 


No. 


n 


L' 


n 


L 


Vl2 J 


1 


3 


1 


3 


1 


1 


2 


3 


1 


3 


1 


1 1 


3 


3 


1 


3 


1 


2 


4 


3 


1 


3 


1 


2 1 


5 


3 


1 


3 


1 


2 2 


6 


3 


3 


3 


3 


3 


7 


3 


3 


3 


3 


1 3 


8 


3 


3 


3 


1 


2 2 


9 


3 


3 


3 


3 


2 2 


10 


3 


3 


3 


3 


2 3 


11 


3 


3 


3 


3 


2 4 


12 


4 





2 








13 


4 





2 





2 1 


14 


4 





2 


2 


2 1 


15 


4 


2 


2 


2 


2 


16 


4 


2 


2 


2 


1 2 


17 


4 


2 


2 





2 1 


18 


4 


2 


2 


2 


2 1 


19 


4 


2 


2 


2 


2 2 


20 


4 


2 


2 


2 


2 3 


21 


4 


4 


2 


2 


2 3 


22 


5 


1 


1 


1 


1 


23 


5 


1 


1 


1 


1 1 


24 


5 


1 


1 


1 


2 




c 



1 

X 


1 

L 


X 


Z X 


26 


5 


1 


1 


1 


2 2 


27 


5 


3 


1 


1 


2 2 


28 


6 














29 


6 











2 1 


30 


6 


2 








2 1 



terms: the tensor-like form of the pseudopotential (Ta- 
bles IVIIl IVIIIl and lIXp has more terms than the central- 
like form (Tables HH HVl and|V|. This means that not 
all of the terms of the tensor-like form arc linearly in- 
dependent from one another, even though they are all 
allowed by the symmetries, and thus some terms can be 
expressed as linear combinations of others, or, equiva- 
lently, some linear combinations of terms are identically 
equal to zero. This fact, can be expressed in the form of 
the following explicit dependencies between the parame- 
ters of the tensor-like pseudopotential. 
For the second-order terms we have. 



/=<ii 

"-"11,21 



2 /=.ii 

^^11,20 



5/^11 
2'-"ll,22' 



whereas the fourth-order dependencies read. 



G 



9 

2 p. 



X-.22 
^22,22 



22,21 = 

1^31 _ _ ^ r'Si 

'-'11,21 ^ /77^11,20 



21C. 



22 

22,23' 



^^VSl 
3(^11,; 



and finally at sixth order we have. 



/=<31 
L'31.21 



/=<33 
'-'33,23 



c; 



42 

22,21 



— r^i I 

^(^31,20 + 

^A42 
g "-^22, 22 ' 



5 /=.31 
3'-'31,22j 

9 /=(33 
^'-'33,24' 

— ^f^C 22^231 



'-'11,21 ^ 



V3 



- ^51 
/77'-'ll,20 



^/=.51 



(16) 

(17a) 
(17b) 

(18a) 
(18b) 
(18c) 
(18d) 



III. RELATIONS BETWEEN THE 
PSEUDOPOTENTIAL AND ENERGY DENSITY 
FUNCTIONAL 



The reader might have noticed that the two forms of 
the pseudopotential do not have the same numbers of 



The EDF related to the pseudopotential is obtained by 
averaging the pseudopotential V over the uncorrelated 
wavefunction (a Slater determinant), that is. 



J 



^ = T / drirar^rj ^^F(risiiir2S2^2,'riSi^ir2S2i2)p(riSiti,r'iSiii)p(r2S2i2,r"2S24)> (19) 
•' S1S2 titj 

where the two-body spin-isospin matrix element of the pseudopotential is defined as 

V(r[s\t\r'2s'2t'2,riSitir2S2t2) = (s'lt'i, ^t^lFlsi^i, 52*2), (20) 

I 

and p(risiti, r^s'^i'j) and p(r2S2t2, '"2*2*2)1 ^re the one- see, e.g., Ref. [il|-) In this lengthy calculation, one 
body densities in spin-isospin channels. (For definitions, must consider as intermediate step the recoupling of the 
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TABLE X: Second-order coupling constants of the isoscalar 
EDF (t — 0) as functions of parameters of the pseudopo- 

tential, expressed by the formula C^^^^j^'^j" — j4(aCoo,oo + 

''C'00,20 + cCoo,22 + '^C'11,00 + 6(^11,20 + /C*!!,!! + 5C'll,22)- 





A 


a 





c 


fj 


e 


f 
/ 


9 


^0000,0 
'-^20,0000 


1 

32 


Q 

—0 


— V 


U 


c 



— Vo 


U 


U 


^0000,0 
^00,2000 


1 

16 


Q 


R 
y 


U 





/q 
— Vo 


U 


U 


^1110,0 
^00,1110 


1 

48 


y 








— \J 6 


Q 


U 


Dv 


^1111,0 
'-'00,1111 


1 

48 


Q 


Oy 'J 




Q 

— 


ov <5 


U 




^1112,0 
00, 1112 


1 

48 




5\/5 


1 




375 





3 


^0000,0 

^11,1111 


3 
4 

















1 





^1101.0 
'-^ 00, 1101 


1 

16 


-3 


-^/3 





-5 


73 








^0011,0 
*-^20,0011 


1 

32 


^/3 


5 





73 


-3 








^0011,0 
^22,0011 


1 

16 








1 











-3 


^0011,0 
'-^ 00, 2011 


1 

16 


-^/3 


-5 





v/3 


-3 








^0011,0 
^00,2211 


1 
8 








-1 











-3 


^1101,0 
^11,0011 


3 
4 

















1 






TABLE XL Same as in Table|X]but for isovector EDF = 1), 
according to the formula C^j'^^Li/ = y4(aCoo,oo + feC'00,20 + 

^^00,22 + c'C'11,00 + 6*^11,20 + /Cll.ll + 5C'll,22)- 





A 


a 


b 


c 


d 


e 


/ 


g 


^0000,1 
"-^20,0000 


1 

32 


73 


-3 





73 


-3 








^0000.1 
^^00,2000 


1 

16 


-73 


3 





73 


-3 








^1110,1 
'^OO.lllO 


1 

48 


3 


73 


-2715 


-3 


-73 





2715 


^1111,1 
'-'DO, 1111 


1 

16 


73 


1 


75 


-73 


-1 





-75 


^1112.1 
^^00,1112 


1 

48 


3v 5 


V 15 


-v3 


— 3v5 


— V 15 





V3 


^0000.1 
"-^11,1111 


-173 

















1 





^1101,1 
^-^00,1101 


1 

16 


73 


-3 





-73 


3 








^0011.1 
^-"20,0011 


1 

32 


3 


73 





3 


73 








^0011.1 
^^22,0011 


-1^73 








1 











1 


^0011,1 
'-'00,2011 


1 

16 


-3 


-73 





3 


73 








^0011.1 
'-'00,2211 


173 








1 











-1 


^1101,1 
'-^ll.OOll 


-173 

















1 






relative-momentum operators, so as to recast the gra- 
dients in such a way that each tensor affects only one 
particle at a time jl^. Such recoupling was performed 
with the aid of symbolic programming, and is not, for 
the sake of brevity, reported in this paper. 

For each term of the pseudopotential (ITJ , we can write 
the result of the averaging in the following way, 

l/~inL' \ \ ^ f~in' L'v' J' ,trpn' L'v' J' .t fOI 'l 

^ hL,vi2S hL,vi2S' 2^ ml.nLvJ -^ml.nLvJ ' V"^^) 

where C^i^nLvj and TI^^^Jlv/ denote, respectively, the 
coupling constants and terms of the EDF according to 
the formalism developed in Ref. Since here we treat 
the isospin degree of freedom explicitly, to the notation 
of Ref. [2] we have added superscripts t, which denote 
the isoscalar {t = 0) and isovector {t — 1) channels. 

Once relations (PT|) are evaluated for each term of the 
pseudopotential, all terms of the N'^LO EDF are gener- 
ated, with the EDF coupling constants C^^^^lf becom- 
ing linear combinations of the pseudopotential strength 
parameters C";^ „. Since the pseudopotentials are 

nL,vi2S 

Galilean- invariant, the obtained EDF coupling constants 
obey the Galilean- invariance constraints [2]. Similarly, 
when parameters of the pseudopotential are restricted to 
obey the gauge-invariance conditions defined in Sec. Ill Gl 
the resulting coupling constants correspond to a gauge- 
invariant EDF. 

The 12 second-order isoscalar (isovector) coupling con- 
stants expressed by the 7 second-order pseudopotential 
parameters are given in Table |X] (Table IXII) . Simi- 
lar expressions relating at fourth (sixth) order 45 (129) 
isoscalar and isovector coupling constants to 15 (26) pseu- 
dopotent ial p arameters, are available in the supplemental 
material 1311. 



A. Inverse relations 

In Section |TT] we noticed the fact that once either the 
Galilean or gauge invariance is imposed, the numbers of 
parameters of the pseudopotential are the same, at each 
order, as the numbers of coupling constants of the EDF 
for each isospin. This situation allows us to obtain the in- 
verse relations, namely, expressions relating the coupling 
constants of the EDF to the parameters of the pseudopo- 
tential. For the case of gauge invariance, at second order 
they are given in Tables IXIII and IXIIIl at fourth order in 
Tab les and |XVl and at sixth order in Tables IXVll 
and IXVIII As sets of independent coupling constants of 
the gauge-invariant EDF we selected the ones used in 
Appendix G of Ref. [2]. Note that in each case, the pa- 
rameters of the pseudopotential can be expressed either 
by the isoscalar or by the isovector coupling constants. 
For the case of Galilean invariance, analogous expressions 
are available in the supplemental material [31| . 

B. Constraints on the Energy Density Functional 

The zero range of the pseudopotential is at the origin of 
the specific constraints induced upon the resulting cou- 
pling constants of the EDF. Indeed, elimination of the 
pseudopotential parameters from pairs of relationships 
defined by Tables IXlIllXllIl IXlVHXVl and IXVlHXVnl 
leaves us with sets of linear equations that the EDF cou- 
pling constants must obey. At second order, that is, for 
the standard Skyrme interaction, this fact is well known 
and allows us to express the time-odd coupling constants 
through the time-even ones, see Ref. Q for the complete 
set of expressions. We do not yet know if the analogous 
property may hold at higher orders, because this fact 
crucially depend on the arbitrary choice of the indepen- 
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TABLE XII: Second-order parameters of the pseudopotential 
as functions of the coupling constants of the isoscalar EDF 
(t = 0) when the gauge invariance is imposed, according to 
4 ('^^""'"^'''0 ,i,r'Ooii,o I ^0011,0 I 

- ^l,"t'20,0000 + "'-'20,0011 '^'-^22,0011 "t" 
/.^OOll.O , _^0011,0 
^00,1101 ^ '^'-'11,0011 ^ /'-'00,2011 



the formula C?'?^ „ 

7^1101,0 I x-,1101,0 

on 1101 I ^^11,0011 



A 


a 


b 


c 


d 


e 


f 9 



^00,00 

^20 
'-'00,20 
^22 
'-'00,22 



_ 2 
3 

2 
3 

-4 



^11 

'-'11,20 ^ 
^11 

'-'11,11 ^ 
^11 

'-'11,22 ~ 



10 


-6 
-2^/3 





6 

273 
10 





5 

73 

3 

73 





-V3 



-3 


73 
5 





TABLE XV: Same as in Table IXlIIl but for the fourth-order 
parameters of the pseudopotential, according to the formula 



^0011,1 I r^UU 
^'-'40,0011 + ./ '-'42 



0011,1 ^ 
0011; 





A 


a 


b 


c 


d 


e / 


^31 
'-'11,00 


2 

15 


-6\/T5 


-18715 


-2175 


40 73 


120 


^31 
'-'11,20 


2 

15 


ISV5 


-1SV5 


-7715 


-120 


4073 


^33 
'-'11,22 


8 [T 
3 Y 5 








-1 





2 


^22 
'-'22,00 


1 

9 


-6^/3 


-1873 


-21 


-8715 


-2475 


^22 
'-'22,20 


1 

9V^ 


18^5 


-1875 


-7715 


120 


-4073 


^22 
'-'22,22 


4 [T 
3 V 3 








1 





2 



TABLE XIII: Same as in Table IXlIl but for the isovector EDF 
{t = 1), according to the formula C^'i l^^^ = ^('i<^20,oooo + 

7^0011.1 I ^0011,1 I 1^1101,1 , ^1101,1 I /-^OOll.l I 

'''-^20,0011 + '-^'-'22,0011 + "'-'00. 1101 ' '''-'11. 0011 + /'-'00,2011 + 
^0011,1 s 
y'-'00,2211 j- 





A 


a 


b 


c 


d 


e 


/ 


9 


ri20 
'-'00,00 


2 
3 


273 


6 





73 





-3 





r^20 

'-'00,20 


2 
3 


6 


-273 





3 





73 





'-'00,22 


4 
3 








-273 











73 


'-'11,00 


2 
3 


-273 


-6 





73 





-3 





^11 

'-'11,20 


2 


-273 


2 





73 





1 





^11 

'-'11,11 


4 

V3 














1 








^11 
'-'11,22 


4 

V3 








2 











1 



TABLE XVI: Same as in Table |XlI] but for the sixth-order 
parameters of the pseudopotential, according to the formula 

^n'L' _ ^/ ^2212,0 ,7^3303,0 , ^0011,0 ,7^0000,0 , 

7iL,jJ12S ^ '^'.'^'-'00,4212 + '"-'00,3303 '^'-'00,6211 + "'-'60,0000 
^0611,0 I r^0011,0 X 
^'-'60,0011 + ./ '-'62,0011-'- 

A a b c d e f 

'^22,00 4i 21 -157105 4273 -20875 -56715 
<^22,2o -|\/^ 3721 -975 1877 -87105-24735 
<^22,22 -^ 1 -2 

^33,00 -j5 7T05 45 6 735 -40721 407 7 
Cii,2o -5721 -975 -3077 87T05 -40735 



dent coupling constants that define the Galilean or gauge 
symmetries. 

In the present paper, we derive the set of constraints on 
the EDF coupling constants that can be obtained by in- 
verting the relations for the isovector coupling constants, 
given in Tables IXllIl IXVl and IXVUI This allows us to 



express, at each order, the isovector coupling constants 
through the isoscalar ones. For the case of gauge invari- 
ance, at second, fourth, and sixth order, such relations 
are listed in Tables IXVIIIi IXIXl and Ell respectively. 
For the case of Galilean invariance, analogous expressions 
are available in the supplemental material p31,] . 



TABLE XIV: Same as in Table |XII] but for the fourth-order 
parameters of the pseudopotential, according to the formula 

^fi'V _ ^/ ^2202,0 ,7^2212,0 I ^0011,0 I 7^0000,0 , 
'-'^ ? _ o — '^l"'-'00,2202 + '"-'00,2212 + '^'-'00,4211 + "'-'40,0000 ' 



nL ,vi2 S 
^^0011,0 



^40,001 



j_ /■^0011,0 \ 
1 + /'-'42,001lJ 



A 


a 


b 


c 


d 


e 


f 



TABLE XVII: Same as in Table IXlIIl but for the sixth-order 
parameters of the pseudopotential, according to the formula 

r<n'L' _ ^2212,1 .7^3303,1 , ^0011,1 , ,^0000,1 , 

'-'.- F c — '^^'''-'00,4212 "T "'-'00,3303 '^'-'00,6211 "'-'60,0000 
.0011,1 



^0011,1 



A 


a 


b 


c 


d 


e 


f 



'-'11,00 


2 

15 


1875 


-1875 


-7715 


-120 


4073 





^53 
'-'11,22 


^31 
'-'11,20 


2 

15 


6715 


-307T5 


-3575 


-4073 


200 





^42 
'-'22,00 


/^33 
'-'11,22 


8 fV 
3 Y 15 








-1 








2 


^42 
'-'22,20 


^22 
'-'22,00 


1 

9 


30 


18 


773 


40 75 


8715 





,'^44 
'-'22,22 


^22 
'-'22,20 


1 

9v^ 


6715 


18715 


2l75 


4073 


120 





^33 
'-'33,00 


^22 
'-'22,22 


-177 








1 








2 


^33 
'-'33,20 



16 
3 

2 
7^/3 
_2^ 
21 
16 

yis 

___2_ 

45 
_1 
15 



1 2 

-21 37T05 -4273 567 5 56715 
21 97T05 4273 16875 -56715 
1 -2 

-3735 1573 -67T05 -4077 -4072T 
377 97l5 672T -24735 87T05 
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TABLE XVIII: Constraints on the EDF that is derived by 
averaging the second-order gauge-invariant pseudopotential, 



expressed by the formula Cj^/^'^^y^ 



1101,0 
00,1101 



+ bC 



0011,0 



+ 



1101,1 

00,1101 

0011.1 

00,2011 

0011,1 

00,2211 



^1101,1 
'-'11,0011 



0000,1 

20,0000 

0011,1 

20,0011 

0011,1 

22,0011 



7^1101,0 I ^0000,0 
^^"^11,0011 + ^'-'20,0000 



_,_ /.^0011,0 I ^0011,0 
/'-'20,0011 S'-'22,0011 



a 


b 


c 


1 











1 

V3 











1 











1 

2V3 


1 

2 





1 

2 


1 

2^/3 











1 

V3 








1 







2 

2 



-2 
2 








1 









1 








4 

V3 





1_ 



TABLE XXI: Second-order coupling constants of the EDF 
as functions of parameters of the pseudopotential when the 
gauge and the spherical symmetries are simultaneously im- 

= A{aCi 



posed, according to the formula Cjjj/^„"j^^j'* 



20 I 
00,00 + 



bCfy 



' cCoo,22 + c'C'11,00 



" £^711,20 



/C'il.22) 



0000.0 

20,0000 

-0000,1 

20,0000 

iOOOO,0 

00,2000 

i0000,l 

00,2000 

1111,0 

00,1111 

1111,1 

00,1111 



A 



J_ 

32 
J_ 
32 
J_ 
16 
J_ 
16 
J_ 
48 
J_ 
16 



a 


b 


c 


d 


e 


/ 


-3 


-^/3 





5 


-73 







-3 





v/3 


-3 





3 


^/3 





5 


-%/3 







3 





73 


-3 





3 


5v/3 




-3 


3^3 


-37T5 


73 


1 


V5 




-1 


-75 



TABLE XIX: Same as in Table lXVIlIl but for the fourth-order 

terms, according to the formula C^^^nhtj^ ~ '^('^C'oo'Sm + 
7^2212,0 , ^0011,0 , 7^0000,0 , ^0011,0 , /.^0011,d \ 

"'-'00,2212 + CLyoo.4211 + '^'-'40,0000 + 6'-' 40, 0011 + ./ '-'42.0011 / ■ 





A 


a 


b 


c 


d 


e 


,/ 


^0000,1 
'-'40,0000 


1 
120 


-67T5 


-18\/T5 


-21^5 


-40^3 








^0011,1 
"-'40,0011 


1 
120 


-18\/5 


18^5 


7715 





-40\/3 





^0011,1 
"-'42,0011 
^2202,1 
'-'00,2202 


1 
V3 
1 
9 




-3\/3 






1 





-47T5 




-I2V5 


1 



^0011,1 
'-'00,4211 
^2212,1 
'-'00,2212 


1 
%/3 
1 
9 








-3^3 


1 





-47T5 




4^5 


4 
14 



to the case of Galilean symmetry are given in the supple- 
mental material [31.] . 

When the gauge symmetry is imposed on the EDF 
and the isospin degree of freedom is taken into account, 
we have 8 independent spherical EDF terms at second 
order, 6 at fourth order, and 6 at sixth order. The 8 
corresponding second-order coupling constants can then 
be expressed by the 7 second-order pseudopotential pa- 
rameters. Similarly, both at fourth and sixth orders, 6 
coupling constants can then be expressed by 6 pseudopo- 
tential parameters. 

As is well known, at second order the isoscalar and 
isovector spin-orbit coupling constants depend both on 
one spin-orbit pseudopotential parameter, namely. 



IV. RELATIONS BETWEEN THE 
PSEUDOPOTENTIAL AND ENERGY DENSITY 
FUNCTIONAL WITH CONSERVED SPHERICAL 
SYMMETRY 



^0000,0 _ '^^11 

*-'ll,llll — ^'-'ll,llJ 

^0000,1 _ _"v/3x-.ii 

'-'11,1111 — ~"-^ii,iii 



(22a) 
(22b) 



In this Section, we assume the spherical, space- 
inversion, and time-reversal symmetries of the EDF, see 
Sec. IV of Ref. [5]. In this way we make our results appli- 
cable to the simplest case of spherical even-even nuclei. 
Below we fully show explicit results for the case of gauge 
symmetry conserved, whereas the full results pertaining 



TABLE XX: Same as in Table IXVIIII but for the sixth-order 
terms, according to the formula Cj^/„"^^j'^ = A{aCQQ^l2i2 + 

7^3303,0 , ^0011,0 I 7^0000,0 , ^0011,0 , r^OOUfi \ 
O'-'nn MMflM -r COqq gjll + '^'-'60,0000 ^'-'60,0011 / '-'6 



-'62,0011 V 



A 


a 


b 


c 


d 


e 


/ 



^0000,1 
'-'60,0000 


1 
840 


21^15 


-45^7 


126x/5 


-280^/3 








^0011,1 
"-'60,0011 


1 
840 


-21^5 


-45\/2T 


-427T5 





-280^3 





^0011.1 
'-'62,0011 


1 








1 








-1 


^3303,1 
"-'00,3303 


1 
9 





-3^3 





-8^7 


"872T 





^0011,1 
'-'00,6211 


1 








1 








-4 


^2212.1 
'-'00,4212 


1 
3 


"73 








8%/T5 


-875 


-24 



which gives one constraint on the spin-orbit coupling con- 
stants, 



^0000.1 _ 1 /-,0000,0 

'-'11,1111 — "-^11,1111- 



(23) 



The remaining 6 spherical EDF coupling constants ex- 
pressed through 6 pseudopotential parameters are given 
in Table IXXH Similar expressions at fourth and sixth or- 
ders are given in Tables|MIIland|20aiIl As in SecHEl 
from these results we can obtain the inverse expressions 
relating the parameters of the pseudopotential to the cou- 
pling constants of the spherical EDF; these are given in 
Tables IXXTVllXXVIl 

At second order, the gauge and Galilean symmetries 
are equivalent to one another When at higher or- 
ders the Galilean invariance is imposed on the spherical 
EDF, we have at fourth (sixth) order 18 (32) independent 
terms, of which 4 (8) are of the spin-orbit character. It 
turns out that, in the same way as for the second or- 
der, the higher-order spin-orbit coupling constants are 
related only to the spin-orbit pseudopotential parame- 
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TABLE XXIL Same as in Table IXXj but for the fourth- 
order coupling constants of the EDF, according to the ior- 
mula C"'/''"^'^}' - ''^"'^31 ' ' -^^^ ' -J-^^^ 

6^22,20 + /C22,22 



AiaCitoo + bCii^20 + cCii,22 + ^^^22, 00 



A 


a 


b 


c 


d 


e 


f 




TABLE XXV: Fourth-order parameters of the pseudopoten- 
tial as functions of the coupling constants of the EDF when 
the gauge and the spherical symmetries are simultaneously 
imposed, according to the formula C" ^ 



^0000,0 I 



66, 



AL,vi2S — "-^40,0000 
0000.1 , ^2202,0 , j^2202,l , ^1111.0 , j.^1111,1 
40,0000 
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75 
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16 
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3 
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64 
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r22 
■22,00 


8\/5 
3 


8 /s" 
~3Y 3 


2 


2 


t22 
22,20 


8 /i" 

3 Y 3 


8\/5 
3 


2 

Vz 


2 


r22 


32 


32 /s" 


4 


4 


■22,22 


3 Y 21 


3 Y 







e 


/ 














4 


4 




721 














-2. 

V 21 





TABLE XXIII: Same as in Table IXXj but for the sixth- 
order coupling constants of the EDF, according to the for- 



mula Cii^'^:ii'f = A{aCr 
eC7,33 . ^^33 



1,22 + ''C'22.00 + cC22,20 + dC22,22 + 



33,00 



0000,0 1 

60,0000 4480 

0000,1 1 

60,0000 4480 

0000,0 1 

00,6000 2240 

0000,1 1 

00,6000 2240 



-21^-7\/T5 
7VT5 -2175 
21V5 7VT5 

-7^/T5 21^ 



d 


e 


f 





75V2I 


-45\/7 





45V7 


-45V2T 





75V2I 


-45\/7 





45V7 


-45V2T 



Coo,3i'ii m -135\/7 21^5 35\/T5 -10573-45^ 135^ 
^ I5V2T -7/15 -7^5 -105 4577 15^21 



r" 



00,3111 800 



TABLE XXVI: Same as in Table IXXVl but for the sixth-order 
parameters of the pseudopotential, according to the formula 



^h' L' 
iiL .V12 S 



00,3111 



_ ^0000,0 
— '^'^60,0000 
I /•^3111.1 
^ 7^00,3111 • 



. ,^0000,1 . ^0000,0 
+ "'-'60,0000 + '^'-'00,6000 



,^0000,1 
' "'-'00,6000 + 



^53 
(-^11,22 



'22,00 



'22,22 



-'33,00 
33 

33,20 



64^ 
9 

-16^5 

64 
3%/3 
16 /T 
3 Y 3 
16^ 
9 



64 Tt~ 
""9"Y 3 

-16\/5 

_ 64 
3 

16v/7 
9 

_ 80 
9 



16^ 
9 

875 



5 
3 

16 
3v^ 



8\/5 

_ 16 
3 

-177 



8^ 
9 



40 
9v^ 




40 
' 21V3 






40 

' 9V2T 




40 
21 






ters. Namely, at fourth order we have 



TABLE XXIV: Second-order parameters of the pseudopoten- 
tial (spin-orbit term not included) as functions of the coupling 
constants of the EDF when the gauge and the spherical sym- 
metries are simultaneously imposed, according to the formula 



^1111,0 
f^^nn 1111 



'20,0000 

/•^1111,1 

" /'-'00,1111 • 



-I- oG. 



0000,1 
20,0000 



0000,0 
00,2000 



, ,^0000,1 
+ "'-'00,2000 + 



(-'( 



20 

00,00 



^22 

'-'00,: 

^11 

'-'11, ( 

^11 

'-'11,: 

^11 

'-'11,: 



-4 ^ 

4_ 

16 

4 

4 
16 



b 


C 




4 


2 


2 


%/3 




73 


-4 


2 

v/3 


2 


16 


4 


4 




Vis 




4 


2 


2 








20 


2 


10 


3 




3 


16 


4 


4 


3^5 




3^5 






12 





4 



^0000,0 
'-'31,1111 



^0000,1 

'-'31,1111 



^0000,0 
'-'11,3111 



^0000,1 _ 
'-'11,3111 ^ 



3 —31 1 1'^ r<'^2 

rj^22,ll; 



/^31 
'-'11,11 



16 
1 

16 



3 --31 1 /^r'22 

g'-'22,lli 



_ " ^31 _ 
-j^gVOOii 11 



^22 



(24a) 
(24b) 
(24c) 
(24d) 



"which gives the following constraints on the spin-orbit 
coupling constants: 



^0000,1 

'-'31,1111 



^0000,1 
'-'11,3111 



1 

2 

7! 



^0000.0 

'-'31,1111 



^0000,0 
'-'31,1111 



2 

1 

73 



^0000.0 

'-'11,3111' 



C^iTsm, (25b) 



(25a) 



and at sixth order we have 



-0000,0 
51,1111 



_3 

64 



" ^51 , 1 /«J^42 ^ ^31 

^•-^11,11 I 7'-"22,ll ~ ^"-"31,11 



3 

64' 



27 / 1 ^33 
■^V 14^33.11, 



(26a) 



c; 



'0000,1 
51,1111 



"-^11,11 32 V c; '-'22,11 "-^31,11 



64 



64 



^0000.0 
^-^11,5111 



^0000,1 
*-^ll,5111 



9 / 3 ^33 
■goV 14^33,11> 

3 --,51 1 /3^42 3 



(26b) 



64 



32 



c^'-^ll,ll on\l ^"^22, 11 64^-^'" 



27 / 1 ^33 
"^V 14^33.11, 



(26c) 



V^,^51 , 3 /^^42 V'-'^3i 
"J7'-'ll,ll + ■^\/ ■^'-'22,11 ~ "Fr'-"31,11 



-0000.0 
31,3111 



64 
_9_ 
80 
21 

160 
27 
200 



64 



^33 
^^'-'33, 11' 

9 

160' 



7^33 
2'-'33,ll' 



/^51 I /^31 

T^'-'ii.ii + T^*"3i,ii 



^0000,1 
•-^31,3111 



— \/3r^i I 



-0000,0 
33,3313 



-0000.1 
33,3313 



9 /21 

~20oVT 

24 V 2 
3 

+ 80 

1 /z 

24 V 6 



^33 

^33,11' 



-C: 



^51 
^11,11 



33 

33,11' 



^51 
^11,11 



1 

24 



1 n 

24 V 6 



^31 
^31,11 



^31 
^31,11 



80 "3341' 
which gives the constraints: 



0000,1 

51,1111 



_ 1 ^0000,0 
^^51,1111 



^0000,1 
•-^11,5111 



^0000,1 
•-^31,3111 



^ ^0000.0 
^^11,5111' 

1 



^ ^0000.0 _ ^ ^0000.0 
_yg^51,llll ^'-'11,5111' 



1 

7!^ 



^ ^0000,0 

/:t^31,3111' 



0000,1 
33,3313 



1 ^0000,0 
^'-'33, 3313- 



(26d) 
(26e) 
(26f) 
(26g) 

(26h) 

(27a) 
(27b) 
(27c) 
(27d) 



If now we consider the Gahlean-invariant and spher- 
ical EDF without spin-orbit terms, we obtain at fourth 
(sixth) order 1 (2) possible constraints among the remain- 
ing 14 (24) coupling constants related to the remaining 
13 (22) parameters of the pseudopotential. These re- 
sults are available in the supplemental material (3l| . Of 
course, such constraints can be imposed in very many 
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■□ Galilean-invariant EDF 
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FIG. 2: (Color online) Number of terms of the spherical EDF 
that is related to a pseudopotential (solid lines). Full squares 
and circles show results for the Galilean and gauge invariance, 
respectively. For reference, dashed lines with open squares 
and circles show the corresponding results for the general 
spherical EDF studied in Ref. 



different ways. We have checked that, in fact, not any 
of the 1 (2) coupling constants of the fourth (sixth) or- 
der spherical EDF can be considered as being dependent 
on all the other coupling constants. In the supplemen- 
tal material we present one example of a possible choice, 
whereby at fourth (sixth) order the coupling constants 

^22 1111 (^00 3111 ^^'^ ^00 3111) selected to be depen- 
dent. A comparison between the numbers of terms of 
the Galilean-invariant and gauge- invariant spherical EDF 
with and without constraints coming from the reference 
to the pseudopotential is plotted in Fig. [5J 



V. CONCLUSIONS 

In summary, in this work we derived the Galilean- 
invariant nuclear N'^LO pseudopotential with derivatives 
up to sixth order and found the corresponding N'^LO 
EDF, which was obtained by calculating the correspond- 
ing HF average energy. Owing to the zero range of the 
pseudopotential, the number of terms thereof is twice 
smaller then that of the most general EDF. We found ex- 
plicit linear relations between the parameters of the pseu- 
dopotential and coupling constant of the EDF. These lin- 
ear relations constitute a set of constraints, which allow 
for expressing one half of the coupling constants through 
the other half. As an example of such constraints, we 
have derived linear relations between the isoscalar and 
isovector coupling constants. The gauge-invariant form 
of the pseudopotential was also derived, and all deriva- 
tions were repeated also for this case. 

We have also analyzed properties of the EDF restricted 
by imposing the spherical, space-inversion, and time- 
reversal symmetries, which are relevant for describing 
spherical nuclei. In this case, by relating the EDF to 
the pseudopotential, at second, fourth, and sixth order 
one reduces the numbers of coupling constants only from 
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TABLE XXVII: Number of terms of different orders in tlie pseudopotential ([2]) and in the EDF up to N'^LO, evaluated for the 
conserved Galilean and gauge symmetries. The last four columns show the number of terms in the EDF evaluated by taking 
into account the additional constraints coming from the relation of the EDF to pseudopotential. 

Pseudopotential EDF 

Not related to pseudopotential Related to pseudopotential 

General Spherical General Spherical 



Order 


Galilean 


Gauge 


Galilean 


Gauge 


Galilean 


Gauge 


Galilean 


Gauge 


Galilean 


Gauge 





2 


2 


4 


4 


2 


2 


2 


2 


2 


2 


2 


7 


7 


14 


14 


8 


8 


7 


7 


7 


7 


4 


15 


6 


30 


12 


18 


6 


15 


6 


15 


6 


6 


26 


6 


52 


12 


32 


6 


26 


6 


26 


6 


N^LO 


50 


21 


100 


42 


60 


22 


50 


21 


50 


21 



8, 18, and 32 to 7, 15, and 26, respectively. Such reduc- 
tion has two origins: (i) at each order 1, 2, or 4 spin- 
orbit isovector and isoscalar coupling constants become 
dependent on one another and (ii) at fourth and sixth 
order one or two non-spin-orbit coupling constants be- 
come linearly dependent on the remaining 13 or 22 ones, 
respectively. Therefore, in spherical magic nuclei one can 
expect relatively small effects related to imposing on the 
EDF the pseudopotential origins, whereas this may have 
much more important consequences in deformed, asym- 
metric, odd, and/or rotating nuclei. We also note that for 
the EDF related to pseudopotential, imposing the spheri- 
cal symmetry does not change the numbers of idependent 
coupling constants as compared to the general case. 

Table IXXVIII gives an overview of the results by show- 
ing the number of terms of pseudopotential and EDF 
with Galilean or gauge symmetries imposed. 
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Appendix A: Time-reversal invariance and 
hermiticity of the pseudopotential 

The pseudopotential studied in this work is a contact 
interaction built with derivative and spin operators. Fur- 
thermore, the choice concerning the formalism is the use 
of the spherical tensors. Under these assumptions, the 
general structure of the pseudopotential is based on the 
following building blocks. 



[K'uK^l]sSv^.s Ji2(rir^;rir2). (Al) 



The final coupling to a scalar ensures that Vq is invari- 
ant under space rotation Moreover, provided that h' ~\-n 
is even, it is also invariant under space-inversion. Now we 
proceed to explore another fundamental symmetry, the 
time-reversal, and later we also require the hermiticity of 
the pseudopotential. 

The time-reversal operator T = —iayK, where K is 
the complex conjugation in space representation, can be 
explicitly applied to the spherical-tensor representations 
of momentum and spin operators, Eqs. (O, ([S]), and ((B]), 
which gives the generic result for spherical tensors. 



(A2) 



where Ta are numerical phase factors. In our case, 
we obtain Tk — —1 for the momentum operator and 
Tcr^, = (—1)^ for the scalar {v = 0) and vector {v = 1) 
spin operators. Moreover, since the Clebsh-Gordan coef- 
ficients are real, rule (|A2p propagates through the angu- 
lar momentum coupling, that is, if phase factors Ta and 
Ta' characterize tensors Ax and Ay , respectively, then 
the coupled tensor, 



^A'V" 



[AxA 



(A3) 



is characterized by the product of phase factors Ta" = 
TaTa' (cf. Appendix B in Ref. (I). Therefore, the cou- 
pled operators appearing in Vq (|Aip are characterized by 
the following values of phase factors. 



Tk^l = (-1)"' 



(A4a) 

(A4b) 
(A4c) 



Finally, because the Dirac delta is real, for Vq we have, 
Ty^ = (_i)"'+»+"i2^ (A5) 

and by taking into account the space-inversion invari- 
ance, it boils down to 



(A6) 
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This justifies the phase factor i^'^'^ in the definition of the 
pseudopotential in Eq. which ensures that for real 
parameters, all terms of the pseudopotential are time- 
even. 

Now we can proceed to calculate the adjoint of the 
operator Vq ()A1|) multiplied by the phase factor derived 
above, that is, 

(A7) 

where we treat the space derivatives of the Dirac delta 
like ordinary numbers and the space variables had to be 
exchanged, r[ o ri and -O- r2. 

Properties of generic spherical tensors under the com- 
plex and Hermitian conjugations are given by the follow- 
ing rules. 



4^ = Ha{-1)^-''A, 



(A9) 



where the phase factors Pa and Ha can be directly de- 
rived from definitions ([3]), (O, and (|6]), that is, Pk = — 1 
and Ha-^ = +1. These rules also propagate through 
the angular momentum coupling, that is, Pa" — PaPa' 
and, for commuting operators, which is the case here, 
Ha" — HaHa' ■ Therefore, we have, 

P\K.,„K'.] =(-l)"+" = +l, (AlO) 



and 



H 



+1. 



(All) 



A 



PAi-l) 



-fix. — I 



(A8) Finally, the adjoint operator of Eq. (jA7l) is given by 



<^i2(rir^;rir2) 



(A12) 



where the last equality results from flipping the order of 
coupling of the operators Kf^,i, and K'-j, which brings 



out the phase factor of (— 1 



(-l)-5. Therefore, 



the time-even tensor i'"^^Vo is not self-adjoint, but we can 
hermitize it by using the expression given in Eq. ([2]). 



40 

00,00 



Appendix B: Relations defining the gauge-invariant 
pseudopotentials 



As discussed in Section III C[ when the gauge invari- 
ance is imposed on the pseudopotential, one obtains a 
specific set of constraints on the parameters and terms 
of the pseudopotential, which result from the condition 
of Eq. dUl). 



^40 
"-^00,20 



^42 
'-^00,22 



^31 
^11,22 



^20 
'-'20,00 



/-.20 



p2 



2v^ 
3 

3 



^22 
"-^22, 00; 



^22 
"-^22, 20' 



X-.22 
_yy'-'22,225 



^33 
"-^11,22' 



\/5^22 
2 "-'22,00' 

V^r^22 
2 "-^22,205 

99 • 



(Bl) 
(B2) 
(B3) 

(B4) 

(B5) 

(B6) 
(B7) 



At fourth order, the gauge symmetry forces seven pa- 
rameters of terms listed in the Table HVl to be specific lin- 
ear combinations of the four independent ones, namely. 



At sixth order, imposing the gauge symmetry forces 16 
terms of the pseudopotential listed in Table Ivl to be spe- 
cific linear combinations of 6 independent ones, namely. 
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00,00 



^60 
"-^00,20 



ri62 

•-^00,22 



^51 
"-^11,00 



^51 
"-^11,20 



-i44 



4^^22.00, 
^ /-.42 

4 ^22,22; 

9 /3^33 
2 V 7 33,00' 

2 V 7 33,20' 



(B8) 
(B9) 

(BIO) 
(Bll) 
(B12) 



^31 

^31,22 



^33 
•-"31, 



22 



9 

To 
9 
5 



TK\ "r'-^ii,22' 



" ^53 

r- "-^11.22' 



33 

33,22 



- ^53 
T7<-'11,22- 



(B21) 
(B22) 

(B23) 



L'11,22 



40 



G 



^40 
^20.20 



C. 



42 

20.22 



,^40 
"-^22, 22 



^42 



/^3i 

'-'31 



31 

31,00 



^•31 
'-'31 



31,20 



9 /_3 

2 V 35 



n\l ^C'll,22' 



7 



4\/5 
7 

4\/5 



^42 
^-^22, 00' 



^42 
^-^22. 20' 



3\/5^44 



-a 



21 

4V5 



22.22' 



,^44 
^-^22, 22' 



3 W — Cn 9 r 



10 



10 



21c 



33 

33,00' 



21C: 



33 

33,20' 



(B13) 
(B14) 
(B15) 

(B16) 
(B17) 

(B18) 
(B19) 
(B20) 



Appendix C: Relations between the central-like and 
tensor-like pseudopotentials 



In the following we present the recoupling formulae 
which connect the two alternative forms of the pseudopo- 
tential of the Eqs. (H]) and (O. We have, 



^£'.5 - -^^^(i-i^..,..)V25TT 



yJ=\L'-Vi\ 

L'+vi 

+ (-i)'^+"'+^\/2jTT 

J=|L'-i.i| 

L+vi 

+ (-i)'^+"'+^'V2jTT 

J=|L-Dl| 

X (l-P*^P'^P^) <5i2(r'ir^;rir2). 



L' vi J 
V2 L S 

L' j] 

V2 L S j 

L vi j] 
V2 L' S\ 

V2 L' S\ 



K' - cr^l)] \k -cr(2) 



n' L' ^1 



nL ^1 



(CI) 



Analogously, the recoupling formula which allows to express the tensor-like pseudopotential through the central-like 
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one reads, 



L' + L 



S=\L'-L\ 



X (-1) 



S+J+vi+L 



L' L S\ 

V2 Vi J \ 



+ (-1) 



J+V2 + L 



L' L S\ 

V2 Vi J \ 



"V2 "vi 



+ (-1) 



S+.J+V1 + L 



, \ L L' S 

V2 Vi J 



+ (-1) 



J + V2+L' 



L L' S\ 



\V2 Vi J 

X (l-P*^P-p-)^i2Kr^;rir2 
I 



V2 ^Vi 



(C2) 



According to the recoupling of the Eq. (jCl|) . we give 
the hst of the relations between the parameters of the 
two forms of the pseudopotential. For the second order 
terms we have, 



^31 
^11,22 



3 2V3 



ri3i I l/=f3i 

'-'11.21 "T ^"-^11,221 



C33 /-^33 
11,22 — ^-^11,227 
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'-'11.22 



/=t20 
'-'00,00' 
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/=.ll 
'-'11,111 
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^-(7,11 



(C3) 
(C4) 
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~C'ii22' (C7) 



(C8) 
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^20 



'-'20,00 — '-^20,00' 



^20 _ /=y20 

'-'20,20 — '-'20,211 



22 _ A=i22 
20,22 — '-^20,21' 



22,00 ~ ^22,02' 



(C18) 
(C19) 
(C20) 
(C21) 



at the fourth order, 



^31 
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10 '''■^^ 2^3 
at the sixth order. 
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y-^62 /-^62 
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/^51 /^oi 

"-^11,00 ~ "-^11,011 



'51 



(C28) 
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1 X. 



C51 ^ /^51 ^ /^51 I /^51 
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Coo /^OO 
31,22 — '-'31,22 



^33 



'-'33,00 — '-'33,03) 



^ n^S I -'^ /^33 



2^-33,22 "^^33,23 y 7^33,24) (C48) 



3^c 



(C46) 



(C47) 



1 



1 



^42 _ ^ /=<42 I /V42 I / ' >=V42 ('r"^S'l 

'-'22,20 — _yg'-'22,21 + _^'-'22,22 i V '-'22,23) l,>-'-JOJ 



7 

15 



/^33 /^33 

'-'33,11 ~ '-'33,13) 



(C49) 



^42 
'-'22.22 



^42 



/=t42 

^-•99 



j^Q '-'22,21 



\/7_^_42 



2V3 



22,22 



(C39) 



^2^,23) (C40) 



^33 
'-'33,22 



^33 /^33 
'-'33.22 ^'-'' 



2 



2V2 



C||,4. (C50) 



33,23 "T 2-v/T4 '^"^ 
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